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Abstract. Open chaotic systems are expected to possess universal transport 
statistics and recently there have been many advances in understanding and 
obtaining expressions for their transport moments. However when tunnel 
barriers are added, which represents the situation in more general experimental 
physical systems much less is known about the behaviour of the moments. By 
incorporating tunnel barriers in the recursive semiclassical diagrammatic approach 
we obtain the moment generating function of the transmission eigenvalues at 
leading and subleading order. For reflection quantities quantum mechanical 
tunneling phases play an essential role and we introduce new structures to 
deal with them. This allows us to obtain the moment generating function of 
the reflection eigenvalues and the Wigner delay times at leading order. Our 
semiclassical results are in complementary regimes to the leading order results 
derived from random matrix theory expanding the range of theoretically known 
moments. As a further application we derive to leading order the density of states 
of Andreev billiards coupled to a superconductor through tunnel barriers. 



PACS numbers: 03.65.Sq, 05.45.Mt, 72.70.+m, 73.23.-b, 74.40.-n, 74.45.+C 
1. Introduction 

Quantum systems that are chaotic in the classical limit exhibit universal behaviour, 
for example for the transport through quantum dots [U |2]. For open systems, the 
transport properties are encoded in the scattering matrix connecting the incoming 
and outgoing wavefunctions in the leads. If we imagine a chaotic cavity attached 
to two scattering leads carrying Ni and iV2 channels respectively, with a total of 
N — Ni + N2, then the scattering matrix separates into four blocks 

siE) = ( ; ) . (1) 

Of particular interest are the transmission eigenvalues and their moments Tr [i^^]" 
given in terms of the transmission subblock of the scattering matrix. The transmission 
eigenvalues relate to electronic transport through the cavity, like the conductance 
which is proportional to their first moment [Sj IH [5] and the power of the shot noise 
which is related to their second. 

To obtain a handle on the transport moments, one can employ the semiclassical 
approximation for the elements of the scattering matrix [6l [7] 

So^iE)^^ E A(^^)e^"^(^\ (2) 
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Figure 1. The trajectory quadruplet which meets at a central encounter in (a) 
contributes to the leading order term of the second moment. The quadruplet 
can be redrawn as the tree in (d) where the paths around the trees recreate 
the trajectories in (a). By moving the encounter to the incoming lead, with 
il = 12, we obtain the possible trajectory configuration in (b), here shown with 
the trajectories tunneling through the barrier. Removing the paths on the left of 
the encounter can be represented by leaving stubs in the tree diagram as in (e). 
Moving the encounter to the outgoing lead provides the structure in (c) or the 
tree diagram in (f). 



in terms of all the classical scattering trajectories 7 that connect the corresponding 
channels. They contribute with their stability amplitude and a phase involving 
their action S-y while /i is the escape rate of the corresponding classical system. 

The first moment of the transmission eigenvalues then follows from a sum over 
pairs of trajectories which start in the same channel in one lead and end in the same 
channel in the other 

Ni N2 

T^N-^EE E A,A;,e^i^^-^^'). (3) 

i=l 0=1 'y(i-S.o) 

7'(i-J-o) 

This fluctuates as the energy is varied and as we are particularly interested in the 
statistics of the transmission eigenvalues, we average over a range of energies. The 
averaging over the phase difference picks out pairs of trajectories with highly similar 
actions and indeed when the trajectories are identical (7 = 7'), one recovers the 
contribution N1N2/N [51 H] which is the leading order term of an expansion in 
of the full result. The next order term was discovered to follow from trajectories that 
come close to themselves in an 'encounter' as they travel through the cavity [TD]. A 
partner trajectory can then be found which is nearly identical, but which traverses 
the encounter differently. Together these pairs provide the weak localisation correction 
—NiN2/N'^ for systems with time reversal symmetry. Higher order corrections came 
from more complicated pairs of trajectories which are nearly identical in long stretches 
called 'links' while differing in small encounter regions and which could be summed 
to give the complete result for the first moment [llj . 

The semiclassical treatment of the first two moments and other correlations 
functions [11] [121 E] led to diagrammatic rules where the contribution of any diagram 
could be read off from its structure. In particular, each link provides a factor of 
A^~^, each encounter a factor of — and each channel in the first or second lead a 
factor of Ni or The main leading order diagram for the second moment of the 
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transmission eigenvalues, which is pictured in figure [Ija) then provides a contribution 
of —N^N2/N^- From this diagram it is further possible to shift the encounter to 
the left until it enters the incoming lead and the incoming channels coincide, as 
in figure [Ub). This configuration provides a contribution of A^iA^I/iV^ while the 
corresponding case of moving the encounter into the outgoing lead in figure W[c) 
provides the contribution NiN2/N'^. These three cases are the only possibilities at 
leading order, giving a total of N^{1 - ^) with ^ = 7Vi7V2/iV^. 

With the diagrammatic rules, the problem of calculating the moments reduces 
to that of finding all the possible diagrams, which was performed for the first two 
moments [11] [121 US] through the connection to correlated periodic orbits which 
contribute to spectral statistics of closed systems [Mj [151 [S] ■ For higher moments 
instead, the leading order diagrams can be represented as trees [T7]. For example, 
the trajectories in figure [Ija) become the tree in figure [Ijd) while moving encounters 
into the leads corresponds to removing alternating links around the encounter node 
as in figures [TJe) and (f). Since trees can be generated recursively by cutting them at 
nodes into sets of smaller trees, all the moments of the transmission eigenvalues |17| 
and the moments of the delay times [18] could be obtained at leading order in inverse 
channel number. This approach could also be adapted to include energy dependence 
and obtain the leading order behaviour of the density of states of Andreev billiards 
[m HO] . Graphical recursions could then be utilised to obtain all moments of these 
various quantities at the next few orders in inverse channel number [21| . 

Exploring the combinatorial interpretation of the semiclassical diagrams, it was 
recently shown that they always combine to give exactly the same results as those 
derived from random matrix theory (RMT) where the scattering matrix is modelled 
as having random elements and belonging to one of the circular ensembles |23| . This 
equivalence between semiclassics and RMT for transport was originally put forward 
in the late 80's [24l [25] and RMT was first used to calculate the conductance and its 
variance |26|, 127) with a diagrammatic approach later providing several leading and 
subleading order results [28]. However, RMT also provides the probability distribution 
of the transmission eigenvalues |23| and their moments are given by integrals related to 
the Selberg integral. This connection more recently allowed the shot noise power [55] 
and then various fourth moments to be calculated [30] and brought a lot of interest to 
the evaluation of these integrals and the corresponding transport moments. Various 
techniques were developed to obtain all the moments of the transmission eigenvalues 
[311 [311 [311 [34] and of the conductance and shot noise [32l|35l|36l [37l [38] as well as of the 
Wigner delay times and time delay [331 13H] ■ Interestingly, the different techniques tend 
to lead to different expressions for the moments which are not so obviously related to 
each other even though they must be equivalent. Similarly, the RMT and semiclassical 
moments must agree [22j but the resulting formulae are different enough to obscure 
the equivalence. However, asymptotic analyses [39] of the particular expressions for 
the moments obtained in [33] have managed to recover the semiclassically calculated 
moment generating functions at the first few orders in inverse channel number [21] . 

The above discussion was for the particular case where the leads couple perfectly 
to the cavity, but for the more general and important case when this coupling becomes 
imperfect much less is known. Imperfect coupling often occurs in physical experiments 
in systems from microwave billiards to quantum dots and to make the theoretical 
treatment above of wider practical use we need to expand the framework to include 
the coupling. For this we add a thin potential wall or tunnel barrier at the end of the 
leads so that any incoming (or outgoing) particle is separated into transmitted and 
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reflected parts, which was originally treated seniiclassically in |40j. Here we show how 
this can be incorporated into the current graphical semiclassical framework [3T]. 

On the RMT side the probability distribution of the transmission eigenvalues 
is currently unknown, with the state of the art being a cavity with one perfect and 
one imperfect lead [41]. The moments are likewise generally unknown apart from at 
leading order when the leads are identical [28]. For the delay times, the probability 
distribution is known j42j but the moments have yet to be evaluated. In general only 
the weak localisation correction to the conductance and the universal conductance 
fluctuations [28j along with the subleading contribution to the shot noise power [IS] 
have been derived diagrammatically. But these have also be obtained semiclassically 
by considering the diagrams explicitly [40| |44| [45] a process we now show how to 
perform implicitly. 

Our paper is organised as follows: In section [5] we explain how tunnel barriers 
lead to modifications in the number of semiclassical diagrams and on the level of their 
individual contributions. In sections |3] and U we derive, to leading order, the moments 
of the transmission and reflection eigenvalues. Section [5] is devoted to the moments of 
the Wigner delay times while in [Appendix A we work out the first subleading order 



terms for the transmission eigenvalues. In Appendix B|we consider Andreev billiards 



with tunnel coupled superconducting leads as an important applications. 
2. Semiclassics with tunnel barriers 

Introducing tunnel barriers leads to two main changes in the semiclassical 
diagrammatic treatment of transport. The first is that the contributions of the 
individual parts of the diagrams changes. The tunnel barriers can be treated 
probabilistically in the semiclassical limit so that trajectories have a certain probability 
to pass through each time they hit the barriers in the leads [ID]. In general this 
probability pi can be different for each channel i. For an ^-encounter involving I 
trajectory stretches of length t which are close together, if the encounter hits channel 
i in the leads, the joint survival probability of all the stretches is (1 — PiY- Over time, 
the survival probability of the encounter is e~'^'* with an escape rate of 

N 

Hl = ^Y.^-{l-Pr)\ (4) 

1=1 

where /i is the escape rate of the same system without tunnel barriers in the lead. The 
diagrammatic rules for the contributions of semiclassical diagrams then become HD] : 



Each link provides the factor y — {^^^iPi^ 



• Each Z-encounter the factor — X^iLi 1^(1^ PiY 

• Each channel sum provides the factor Y^^liP'ii^ — PiY oi" X^i^Wi+i Pi (1 ~ Pi)'' 
where j counts the number of trajectory pairs passing through the same channel 
and k the number of pairs reflected. 

For example for the moments of the transmission eigenvalues, a diagonal pair 
starting in a channel in the incoming lead and travelling together to a channel in the 
outgoing lead give the leading order term of the flrst moment 

(Tr[ttt])^ ^-^^lr'-"-+^^'' +Q(l). (5) 
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Figure 2. With tunnel barriers, trajeetory stretches can additionally be reflected 
at the tunnel barriers. When the encounter from figure [Ja) now moves into the 
lead, one trajectory stretch can tunnel through while the other is reflected as in 
(a)-(c). In the tree representation there are four possibilities of removing one link 
(represented as a stub) and reflecting the other (represented by the perpendicular 
bar) as in diagrams (d)-(f). 



With the same tunnehng probabihty p in each channel, this simphfies to pNS^ with 
^ = NiN2/N'^ while the standard second moment diagram of figure [TJa) provides 

- (1 - (1 - pf) N^^j^ = p{p - 2)Ne. (6) 
Moving the encounter fully into the leads as in figures [Ijb) and (c) gives an additional 

The second main change with tunnel barriers is that many additional diagrams 
become possible [301 HI]- Encounters can now partially enter the leads as trajectories 
can be reflected back into the cavity and are no longer forced to leave the system. 
Staying with the simpler example of the moments of the transmission eigenvalues, 
from the diagram in figure [Ua), we could move the encounter into either lead, let 
one trajectory pair stretch pass through the tunnel barrier and one be reflected as in 
flguresma)-(c). From the channel sum of the encounter touching the lead, we obtain 
a factor of p{l — p)Ni or p{l — p)N2 and an additional contribution of 

Ml - p)^^'^ + Ml - p)^.^ - Ml - P)Ne, (8) 

where the factor 2 comes from the two possibilities of which trajectory to let pass 
through the barrier. The total result for the leading order term of the second moment 
of the transmission eigenvalues becomes 

Tr [th] = pN^ip + 2^ - 3pO + 0(1). (9) 

Our aim now is to systematically generate these additional possible diagrams, 
with their semiclassical contributions, for higher moments. 
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3. Moments of the transmission eigenvalues 

As the moments of the transmission eigenvalues involve trajectories starting in one 
lead and ending in the other, the possibilities for moving encounters partially into 
the leads are somewhat limited. They can therefore be treated relatively easily by 
extending the treatment without tunnel barriers and here we mainly highlight the 
changes needed to incorporate them. At leading order in inverse channel number, the 
underlying semiclassical diagrams can be redrawn as trees |17) as going from the top 
lines in figure [1] and figure [5] to the bottom lines. This framework has been further 
developed [TSl [H] and detailed in [5^ . Graphical recursions can be used to go beyond 
leading order and we build on that formalism here using generating functions as defined 
in |21| . In particular we will obtain an expansion for the moment generating function 

oo 

T{s) = 5" (Tr [i^i]") = iVTo + Ti + . . . (10) 

n=l 



3.1. Tree recursions 

We let (j) and (j) be generating functions which count all trees (including their 
semiclassical contributions) starting inside the system from an incoming or outgoing 
direction respectively and not rooted in any channel. We count the trees by their order, 
or their number of leaves or channels, through powers of the variable r. The trees start 
with a link which connects to an encounter of arbitrary size /, while alternating below 
it are I further trees of the same type as the starting tree and I — 1 trees of the 
opposite type. For example, for (p we include the tree in figure [l]^d) with the channel 
at the top removed. Without tunnel barriers we could further remove all the I trees 
of (j) type below the encounter and move this encounter into the outgoing lead, as in 
figure [TJf). Since the start point of (j) is inside the system, the diagram in figure [TJe) 
is not included, though that diagram with the channel oi removed is included in 0. 

With tunnel barriers the new possibilities of the encounter are illustrated for the 
second moment in figure [51 In general, touching the outgoing lead involves removing 
any k of the I trees and allowing the other Z — A; to be reflected back into the cavity, 
as in figure [2jf). The case fc = just means that the entire encounter is reflected 
at the tunnel barrier and is already included in the semiclassical contribution of the 
encounter. The encounter could also move to the incoming lead with the first link 
from the starting point to the encounter reflected and k of the I — 1 trees of type 
(j) removed to tunnel straight into the lead instead. This possibility is illustrated in 
figure [2Kd), while a diagram like that in figure [2Ke) is not included in cf). Again the 
case fc = is already included elsewhere and we obtain the recursion relation 

N oo N 

i=Ni + l 1=2 1=1 



N I , , , 



1=2 i=Ni + l k=l ^ ^ 

+yEE^'E( ^ (11) 

1=2 i=l k=l ^ 

The first line starts with the contribution of the smallest tree which is just a link that 
tunnels through into the outgoing lead and adds all the trees whose top encounter 
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does not move into the lead. The remaining two hnes count the possibihties that 
the encounter moves into the outgoing or incoming lead and k of the stretches of 
the encounter tunnel straight through the lead (with the remaining stretches being 
reflected back into the cavity). Since the powers of r count the order of the diagram, a 
factor r is included with each tree which is removed along with a corresponding factor 
of p to record that those links have tunnelled through the barrier. 

In the recursion relation we can perform the sums over k. Since the (1 — p)' in 
the top line of (|lip corresponds to the fc = terms of both sums, we can simplify to 

N N oo N 

i=l i=Ni+l 1=2 1=2 i=Ni+l 

+ EE'^'(i-p^)(''p^ + '^(i-p^)) ' (12) 

1=2 i=l 

where we further divided both sides by y. Neatly, the first two terms can be rearranged 
into the I = 1 terms of the three sums over Z, so that the recursion reduces to 

_JH_ ^ rpi + - Pt) ^ y> (/)(1-k) ^^^^ 

1-00 i=Ni+i 1 ~ ^Pi^ ~ 00(1 ~ P^^ i=\ 1 ~ '^^''0 ~ 00(1 ~ Pi) 

For <\) we likewise have 

Af0 rp,+0(l-K) ^ ^ 0(1 -P») 



> Pi -t- VK'^ - Pi) ^ H'K'^ - Pt) ^^^^ 

1-00 1 - rpi4> - 00(1 - Pi) i=Ni+i 1 ~ ''P'0 ~ 00(1 ~ Pi) 



3.2. Leading order moments 

To obtain the leading order diagrams, we now root the top of in an incoming channel. 
This also allows us to remove the top link and tunnel straight into an encounter in 
the incoming lead, so we can now add diagrams like figure [Ije) and figure [21[e). Of 
the remaining / — 1 alternating trees emanating from the encounter, any k can also 
move straight into the lead. The generating function Tq, which counts all the leading 
order diagrams with their contributions and divided by N , is then given by 

4=1 i=l 1=2 k=0 ^ ^ 

Performing the sums gives 

NTo ^ y . (16) 

^ l-rp,0-00(l-p,) 



3.3. Fixed tunneling probabilities in each lead 

With the formulae in ([T3| , ([T4| and ([T6| we have a formal generating function for the 
leading order terms of all moments and we can expand 0, and then Tq in powers of r. 
However, with the sums over the different tunneling probabilities in each channel, it is 
difficult to manipulate the expressions further. To proceed we can make the simplifying 
assumption that the channels in each lead have the same tunneling probability of pi 
for the incoming lead and p2 for the outgoing lead. Equations (fT3l) and (fT4| lead to 
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quartic equations for (j) and (f) which in turn lead to the following quartic equation for 
To: 

{s - 1) {s^pIpI + 2spiP2{2 -pi- P2) + (pi - P2f) To 

+2.s(.s - l)piP2 ispiP2 + 2-P1-P2) 

+s [s{s - 1)pIpI+pI{p2 - l)Ci +pI{Pi - 1)C2 

+e {2sp,p2{2-pi -P2) + s{2s - 1)pIpI + [pi ^P2f)] 
WpIpI^{2s - l)To + s^plpie = 0, (17) 
where Ci = Ni/N, (2 = N2/N and = s is the generating variable for the moments 
since the n-th moment involves trees with 2n leaves. Expanding Tq in powers of s we 
need to choose the correct value for the first moment, which is PiP2^/(piCi +P2C2) 
from ([5]). 

3.4- One perfect lead 

Equation (|17p allows us for example to consider the situation where only one of the 
leads has a tunnel barrier by setting pi or p2 to 1, as in the recently considered 
situation where the probability distribution of the reflection eigenvalues was obtained 
from RMT for systems without time reversal symmetry |41j although with different 
tunneling probabilities in the channels in the remaining lead. With equal probabilities 
instead and p2 = 1, we actually obtain a cubic equation semiclassically 

{s - 1) (1 - 1)) r3 + sis - + C2)r2 

+s (C + - 1) + C2(spi - 1)) To + s^pi^C2 - 0, (18) 

which can also be obtained from (flT)) through using ^ = C1C2 and Ci = 1 ~ (2- The 
generating function in (fT8|) should then match the leading order moments found by 
integrating the corresponding transmission probability distribution to that in |41] . 

3.5. Equal tunneling probabilities 

To proceed further we instead set all of the tunneling probabilities equal to p — pi — 
P2. The quartic equation for now has the expansion 

= C2r + e(p + C2(l-2p))r3 + 

+ ^{p^+ K2(l - 2p) + K(3 - 5p) + 2C2e[l - 5p(l ~ p)]) r^... (19) 

while for (p we swap (2 and Ci . The second term in (fT9)) corresponds to the sum of the 
diagrams in figures [D^d) and (f) and figures [SJd) and (f). The quartic for To reduces 
to 

(,s - 1) {sp^ + 4(1 - p)) T^ + 2{s - 1) (sp^ + 2(1 - p)) T^ 

+ [s{s ~ l)p2 + p - 1 + ,se (4(1 - p) + (2s - l)p2)] Tl 

+s(2s - \)v^iTo + sVe^ = 0, (20) 

and again for the expansion in powers of s we need to pick the correct value of for 
the first moment (since we divided To by N^. The first few terms in the expansion are 

To-Ks+pe(p + 2e-3K)s' 

+ K {V^ + 2p(3 - 4p)e + (6 - 21p + 17p2)e^) + . . . (21) 

with the second term corresponding to the second moment calculated explicitly in 
section [2j 
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3.6. Equal leads 

As a further simplification, we can also have an equal number of channels in each lead 
so that Ci = C2 = 1/2. Due to the symmetry we have cj) = (p, so reduces to 

20 rp + 20(l-p) ^2 n ^99^ 

T3^-l_,p^_^2(i_p)' r0-20 + r = O, (22) 

where (j) is actually given by the same quadratic as when there are no tunnel barriers 
in the leads. The moment generating function then satisfies the quadratic equation 

4(s - 1) (s/ + 4(1 - p)) T^ + 4{s - f ).s/To + sV = 0, (23) 

which is also a factor of ([20]) when Ci = C2 = 1/2 or ^ = 1/4. The moment generating 
function can then be given explicitly as 

^ ^ (p(l - g) + (P - 2)VT^s) 
° 2(s-l)(V + 4(l-p)) ' ^ ' 

where we chose the solution of (|23p which matches the first moment p/4. 

3.7. Comparison with RMT 

The RMT result for the leading order probability distribution of the transmission 
eigenvalues was calculated in [28 . However, a final result could be obtained if they 
assumed the two leads were identical. The channels could have different individual 
tunneling probabilities though, as long as the probabilities are matched in the other 
lead. The result was 

Nl fry N 

P{Z) = y ^'^^ =. (25) 

^7r(p2_4p,;Z + 4Z)V^(l-Z) 

Semiclassically instead it is simple to have different sized leads, but with a constant 
tunneling probability in each. To compare the different results we can look at 
the common results for the simplest case of identical leads with a single tunneling 
probability. From the moment generating function in (j24p . including the 0th moment 
as 1, we can perform the Hilbert transformation to get the probability density 

" Jo N 1-sZ ' N 2^(p2_4pZ + 4Z)^Z(l-Z) ^ ^ 

which is exactly (I25p with identical tunneling probabilities. That the distribution in 
(|25p for identical leads is a sum over the channels means that the moment generating 
function would likewise be a sum over terms like in ()24|) again with different pi. 
Semiclassically, with identical leads we have (j) — (j) through symmetry, but this 
result still cannot easily be seen from (fT3|) and (fT6| . Instead we have access to the 
complementary regime of different leads with equal tunneling probabilities. 

We can proceed to higher orders in the inverse number of channels (see jlOl ) 
by similarly modifying the approach 12 Ij without tunnel barriers. For example, we 
present the calculation of Ti in [Appendix A[ However, the more interesting case occurs 
when we consider refiection quantities. 




Figure 3. For reflection quantities incoming and outgoing channels can coincide 
so that with tunnel barriers new possibilities arise for encounter stretches to tunnel 
into the lead. Allowing the stretch to channel 02 from figure [2^) to also tunnel 
into the lead so that channels 01 = 22 give the possible diagram in (a). In 
the graphical representation in (d) this corresponds to removing adjacent links 
around the encounter node, which can be performed in 4 ways. We can remove 
an additional link as in (b) and (e), if 02 = 01 = 12 where now a pair of trajectories 
are directly refiected at the tunnel barrier, and finally have the situation in (c) 
and (f) where none of the trajectories ever enter the system. 

4. Moments of the reflection eigenvalues 

When we consider the moments of the reflection eigenvalues, or their generating 
function 

00 

R{s) = ^ s" ^Tr [rtr]") ^NRo + Ri + ... (27) 

we now have a semiclassical approximation where the trajectories all start and end in 
the same lead. With tunnel barriers, this allows even more diagrammatic possibilities 
when moving encounters into the lead, and we may also have trajectories that never 
enter the system and are reflected instead directly at the tunnel barrier [iOl HI] ■ For 
example for the second moment, along with the diagrams in figures [T] and the 
diagrams in figure [3] are also now possible. 

Along with the additional possibilities, the tunneling phases now become 
particularly important [40j . Looking in detail at the tunnel barrier drawn in 
figure a), where the incoming channel 12 is identical to the outgoing channel Oi, both 
solid trajectories tunnel through while both dashed, complex conjugated, trajectories 
are refiected, one on each side of the barrier. If p is the refiection amplitude of the 
barrier and r the transmission amplitude then the four trajectories of the encounter at 
the barrier give the factor (t/o*)^ which in the semiclassical limit is equal to — p(l ~p) 
with an additional minus sign arising from the quantum mechanical tunneling |40) . 
With equal tunneling probabilities in every channel, the diagrams in figures [Hta) or 
(d) then give a total semiclassical contribution of 

- 4p(l - P)Ni^ = -4p(l - p)NCl (28) 

In figures[3Jb) or (e), the tunneling and reflecting trajectories at the tunnel barrier 
pair a trajectory stretch with a complex conjugated one giving the standard factor 
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of p{l — p) while in figures [3{c) or (f) none of the trajectories ever enter the system. 
Combined they give 

M^'-p)NCi+{l-p)'NC,, (29) 

while the diagrams in figures [T] and [2] provide the contribution 

p{p - 2)NCt + 2p^NCf + 4p(l - p)NCt, (30) 

since now all the channels are in the first lead. We can then, using Ci(l — Ci) = ^, 
write the leading order term of the second moment of the reflection eigenvalues as 

(ti- [rV]^^ = A^Ci +pN^{p-2 + 2^~ 3^) + 0(1). (31) 

Through the unitarity condition, we have 

r'fr + tH = lN„ Tr [rV]^ = TVi - 2Tr [i^t] +Tr (32) 

which is satisfied by the leading order second moments in © and (PT|) since the leading 
order first moment of the transmission eigenvalues from ([5]) is pN£^. 

Auxiliary trees 

For the moments of the reflection eigenvalues, we first let / and / be generating 
functions which count all trees including their semiclassical contributions which start 
inside the system from an incoming or outgoing direction respectively (without being 
rooted in a channel). Due to the symmetry, we actually have / = / but we will treat 
them separately for now. The / trees, which now include diagrams like figures ISJa) 
and (b) with the channel ii removed, initially start from a link connecting to an l- 
encountcr which is followed by I subtrees of type / and l—l trees of type /. In general, 
and as illustrated in figure [31 we can now allow any of the links around an encounter 
node to tunnel straight into the lead and be removed in the graphical representation. 

Counting the possibilities recursively would be straightforward, but for the 
quantum mechanical tunneling phases which actually depend on how many adjacent 
links are removed together. To keep track of these phases we introduce auxiliary 
generating functions wi^a and wi_a which count the contributions below the encounter 
(and not the top link of /) and where the encounter is in a particular channel in the 
lead with tunneling probability p. The subscript a = 1 represents that the last of the 
subtrees tunnels straight into the lead and is removed from the diagram and a = 
represents that it is reflected. In w we ensure that the first subtree tunnels into the 
lead while in w we ensure it is reflected. The smallest encounter is I ~ 2 for which the 
central tree can either tunnel or be reflected so, as in figures |3I a) and (c), we have 

W2,i = Tp*pT*r^ + {TT*)\^f, W2,o = {rp*fr^f + rT*pp*rff, (33) 

and 

W2,i = {pT*fr'^f + pp*TT*rff, W2,o = pT*Tp*rf + {pp*ff^f, (34) 

where for simplicity we have included an encounter entirely reflected at the lead in W2.0 
which needs to be remembered later. With the symmetry we also have w*q = lii^i and 
we can simplify by setting tt* = p, pp* = {1 — p) and {tp*Y = {pT*)'^ = ^^(1 ^p)- 
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(a) I ; (b) I I ; (c) I 




Figure 4. For the auxiliary generating functions w where the first link after 
the encounter node tunnels directly into the lead, we need to track exactly which 
trajectories tunnel through or are reflected at the barrier. For wi i the last link 
tunnels so that, depending on whether the middle link tunnels or is reflected, 
there are two possible configurations in Jii2,i as shown in (a). For wi q the last 
link is refiected so that again there are two configurations in 102,0 shown in (c). 
Starting from all the configurations in wi i and ui; g we can add an additional 
two links on the right of the encounter. If the last tunnels into the lead, the two 
possibilities for the other generate the configurations in ui;_|_i 1 as depicted in (b) 
while if the last link is reflected we obtain wi^_i q and the recursions in Ij35ll . 



Since we know the behaviour of the subtrees at the edges, we can recursively 
generate w/+i,a from wi,a by adding two more subtrees on the right and allowing both 
possibilities for the first. Tracking the trajectories as in figures |4l[b) we obtain 



Wl + 1,0 = Wlfi 



* 2 I ^ r 

pp r + TT rj 



TT*rf + pp*ff 



wi.o 



p*Tp* 



TT*rf 



rf + PP*ff 



(35) 



where the terms with fractions involve replacing a reflecting trajectory by a 
transmitting one or vice versa, and have the values t* pr* / p* = —p and p*Tp* /t* — 
{p — 1). Because this recursion does not depend on whether the first subtree tunnels 
or is reflected, we have the same equations for w. We wish to sum over encounters 
of all sizes, so if we set Wa — X]i^2 ^^'^ likewise for w we obtain the following 
coupled equations 



Wi — p{l — p)r ~ p r f 



Wl 



(1 -p)r^ +prf 



Wopr 



f-r 



Wo - p{l - p)rf{f - r) =wo prf + (l - p)ff +wi{l-p)f f-r 

We obtain a similar equation for w and the solutions, with f — J , are 

pr2[r(l - p- P) -pf] 



(36) 



Wl 



Wo = Wl 
Wo 



1 — (1 — p)r'^ — 2prf — {1 — p — r'^'^p ' 
prf{p- l)(r - /) 



1 - (1 - p)r2 - 2prf - (1 - p - r2)/2 ' 
fHl-p)[f{l-p-r^)+pr] 
1 — (1 — p)r2 — 2prf — (1 — p — r2)/2 



(37) 
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4.2. Tree recursions 

To obtain the generating function / we simply add a link to these contributions and 
sum over all channels in the lead 

J. Ni 00 N 00 N Ni 

y 1=1 i=2 i=l 1=2 i=Ni+l i=l 

where the sum over the (1 — piY terms in the first lead are already included in wq. 
With equal tunneling probabilities in all the channels, this is 

/ /(1-P)[l-Cl/^(1-P)] ^ . . ^ ^ „ ^ . . 

l_ p = 1 _ p^l-p) \-Ci{rp + wi+2wa + wo), (39) 

which leads to a quartic for / whose expansion is 

/ = Cir + e(l-p-Ci(l-2p))r3+ (40) 

+ e ((1 - Pf - Ci(l - P)(l - 2p) - C(l - P)(2 - 5p) + 2Cie[l - 5p(l - p)]) r^ . . 

The second term now corresponds to the sum of the diagrams in figures [TJd) and (f) 
and figures[2fd) and (f) along with two new possibilities like the diagram in figure[3l[d) 
and one diagram like figure [3l[e). 

4.3. Leading order moments 

To move from the generating function / to the leading order moment generating 
function for the refiection eigenvalues Rq we need to add a channel to the top of the 
diagrams in / and also allow the top link in / to tunnel straight into the lead. Since 
we know whether both the outside subtrees of w and w tunnel or are refiected we 
can simply replace their top lead by a tunneling one and replace the corresponding 
transmitting and reflecting eigenvalues of the trajectories 

^1 Ni , ^ ^ ^ * s Afi 

NRo^y^rpJ + y2r{^w^ + ^wo + ^m + ^wo)+y2r'{l~p^), (41) 

7^1 TP TP* pp* J ^ 

where the last term is the contribution to the first moment from a trajectory pair that 
never enters the system. We have pr* /{rp*) = — 1 so that in total we obtain 

NE - V rC^[r{l-p,-p)+p,f] 

° f^^l-{l-p^y-2p,rf-{l-p,-r^)p- ^ ' 

With equal tunneling probabilities, we again obtain a quartic for Rq which can 
be simplified to 



[4(s -l){l-p) + sp"] (kt + Rl + sikt) + sp'Rl 

+ [{s+p- 1)(1 - s + sp) + sp^(\ Rl + s(l + s)p2^7?o + sVC^ = 0, (43) 
where Rq = (1 — s)i?o ~ Ci* removes the C, dependence. The expansion provides 
i?o = Ci(-s + + s^)-p^s+p^ (p _ 2 + 2C - 3K) 

- K (p2 - 3p + 3 - (6 - 15p + 8p2)^ + (6 - 21p + l7p^)C^) + . . . (44) 

The coefficient of is the same as (|5T|) and unitarity can be checked against (PT|) . 
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With equal leads (i ~ 1/2, we end up with the same quadratic rf^ — 2f + r — 
that we had for the transmission, while the moment generating function is 

„ _ s{2-p){{2-p){1-s)+pVT^s) 

° 2(s-l)(4(s-l)(l-p) + sp2) ■ 
Transforming to the probability distribution of the reflection eigenvalues, we have the 
same density as in (|26p but with Z replaced by 1 — Z which is the same mapping as 
between the transmission and reflection eigenvalues. 

4-. 4- One perfect lead 

Since the start and end channels are all in the same lead, is it straightforward to make 
the other lead transparent. Keeping pi — p in the first lead and setting pi = 1 in the 
second means that we remove the sum over the channels in the second lead in (|38| 
and set 1/y = 7V(pCi + (2)- This changes ^ to 

J— =Ci/(l-p) + Ci('^P + w'i +2u.o + wo), (46) 

while remains the same. This leads to a cubic for / and Rq, where again using 
Rq — (1 — s)i?o — Ci^s leads to the simpler form 

(s+p-l)i?3 + sp(l + C2)i?o 

+s[{l-s + sp)C2 + is+p- m Ro + s^pC2^ = 0. (47) 

Of course we could instead allow the first lead to be transparent and set Pi = I 
there and pi — p in the second lead. In this case the encounters can no longer partially 
enter the lead so the number of possible diagrams is drastically reduced to just those 
where the encounter enters the lead fully. We actually then have almost exactly the 
same recursions as when there are no tunnel barriers in either lead, but just with the 
minor corrections to the survival probabilities of the encounters and links. As in |21j 
we obtain 

« 00 00 00 

^ = <i - E + <i E + ^2 E(i - p)'/""'' (48) 
1=2 1=2 1=2 

where the last term is due to the change in the survival probability of the encounters 
while the links provide \/{Ny) = +pC2- Simplifying we get 

/ - <i I (i-Ci)(i-P)/ .49^ 

I-P l-rf^ l-il-p)P ' ^ ' 

where the last term is the correction due to the tunnel barrier. The moment generating 
function is still given by 

Ro - f^:, (50) 
1 - rf 

which leads directly to the cubic 

(s - l)(s + p - l)Rl + s [Ci(3s + 2p-3)-p] Rl 

+<i [Ci (3s + p - 1) - p] Ro + - 0. (51) 

Shifting the generating function as before, Ro — (1 — s)i?o — Ci*, we then obtain exactly 
(|47|) but with C2 replaced by Ci, which leaves ^ unchanged. Swapping (i and C2 just 
means we are considering the moments of the reflection eigenvalues of the second lead 

00 

R'{s) = E Tr [r'V'j " = NR'^ + R[ + . . . (52) 

n=l 
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while the unitarity condition r'r + tt"^ = I^^ ensures that 

R{s) - = R [s) - , 53) 

1 — s 1 — s 

so this much simpler treatment provides the same generating function (|T7|) as the 
full auxiliary tree combinatorics when one lead is transparent. This leading order 
generating function should also arise as the first term of an asymptotic expansion of 
the recently derived RMT probability distribution |41j when the remaining tunneling 
probabilities are set equal. 



5. Moments of the Wigner delay times 



The tree recursions for the moments of the reflection eigenvalues can easily be modified 
to obtain energy dependent generating functions like 



C(e,n) 



1 

N 



■Tr 



sUe 



S\E 



(54) 



which are related to other physical observables like the density of states of Andreev 
billiards and the moments of the Wigner delay times. The generating function GiyS) = 
s"'C(e, n) can be obtained by simply considering the reflection eigenvalues with 
a single lead Ci — 1 a-nd including the energy difference. This changes the encounter 
and link contributions to 



-Nil-il-pY-la), y-'=Nip-a), 
respectively with a = ie. The tree generating function becomes 

/ ... ^ , a/ 



(55) 



f{l-p) + 



+ rp + wi + 2wo + Wo , 



while vj and w remain unchanged. The terms can be combined and simplified to 

/(l-/2-a)_ pr+(l-p-r2)/ 



(56) 



(57) 

1 so that we 



[1 - /2]2 1 - (1 - p)r^ - 2prf - (1 - p - r2)/2 

The leading order generating function Gq is still given by (|42|) with Ci = 
obtain a quartic equation for / and Gq. 

We show how to use the energy dependence to obtain the leading order 
contribution to the density of states of Andreev billiards in [Appendix B| and 
concentrate here on the moments of the Wigner delay times. The delay times are 
the eigenvalues of the Wigner-Smith matrix [46l |47] 



(58) 



and are a measure of the time spent inside the scattering cavity. Their moments can 
be obtained [El [21] through the correlation functions 

'-^s^(e-'J^)s(e+'J^]-i 



Die, n) = — Tr 



(59) 



by differentiating 



1 , 
"n! de" ^ ^ 



(60) 
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If we denote the moment generating function by 

oo 

M(s) = 5]M'V'(Tr[Qn, (61) 

n=l 

then we can obtain an expansion for the leading order lower moments by expanding 
(|59|) binomially and simply substituting the correlation functions C(e, n) 



X 2 2 6 , 2(p+10) 4 10(2p+7) . 

p p^ p-^ p^ 

5.1. Tree recursions 

However to obtain the full moment generating function we need to account for the 
identity matrix in the brackets in ((59|) that creates the difference from (f54|) . As was the 
case without tunnel barriers |18[ 121] , this matrix can be thought of as coming from 
diagonal trajectory pairs which just travel directly from the incoming to outgoing 
channels with no energy difference. With tunnel barriers, such diagonal pairs can be 
formed whenever an encounter goes into the lead and an / subtree remains, as long 
as the other trajectories either side end in the lead so that we have stubs either side. 
The subtree contribution already includes a diagonal pair and to subtract the identity 
matrix we just subtract the contribution of a diagonal pair with no energy difference 
(a = 0). For every / subtree surrounded by two stubs we then replace its contribution 
by / — r /p. This breaks the symmetry between / and / and we look at their two 
generating functions separately. 

The auxiliary tree recursions for /, make it easy to see when a / subtree is added 
with a stub either side. Modifying the recursions appropriately we find 

wi-p^r'^{f -r) =pr{f -r){wi+wo), 

Wo - p{l - p)rf{f - r) ^ wq prf + {1 - p)ff + wi{l - p)f{f ~ r), (63) 
and similar equations for w. This leads to the equation for / 
/(I -//-a) pr+{l-p-pr^)f 



[1 - ff? l+pr^- pr{f + f)-{l-p- pr^)ff - pr^f 



(64) 



For / however, because we add two new subtrees on the right in the auxiliary 
tree recursions, it is not so straightforward to see when a / subtree is surrounded by 
two stubs. Instead we can simply add either a / subtree or a stub to either side of the 
auxiliary trees for /. Doing so we find 

/(!-//-«) ^ pril+r'^) + il-p-pr^)f ^^^^ 

[1 - ffV 1 +_pr2 -pr{f + /) - (I - p - pr^)ff - pr^f' 

which, by multiplying (IMl) by / and (|55|) by /, leads to the simple relation / = {l+r'^)f 
and a quartic equation for / or /. 

5.2. Leading order moments 

To obtain the generating function 

OG 

L(s) = ^s"i?(e,n) (66) 

n=l 
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we then either root / in an incoming channel or tunnel straight into the auxiliary trees 
to obtain 

— = pf - wi ~ Wo - wi + Wo - pr, (67) 

r 

where we included the correction ^rwi/p from the diagonal pair above the auxiliary 
wi trees which have stubs on either side. The last term stems from a pair of trajectories 
that never enter the system. Here the w terms correspond to the auxiliary trees from 
the recursions for / as in (|63p so that finally we obtain 

L„ = VrU-r-2r^f + rff) ^^^^ 

1 + pr^ - pr{f + /) - (1 - p - pr^)// - pr'^f 

and the quartic 

p^{l + s)Ll + 2aps{l + s)(2 - p)Ll 

+ s [{ap — a + p)(a - p) + ap{ap — 2p + 4)s + a^p^s^] Lq 

+ a^p^s^{l + 2s)Lq + aVs^ = 0. (69) 

For the nth moment of the delay times we then differentiate n times (and divide 
by n\) which can be achieved by simply transforming s — > s/a. Setting the energy 
difference a to then provides the generating function 

pMq^ + (4s - 2ps - p)Mo^ + ps^Mn + ps^ = 0. (70) 

Expanding Mq gives the terms in ((62|) directly when we pick the solution whose 
low moments agree with the semiclassical diagrams. The probability distribution of 
the delay times has previously been obtained from RMT including the leading order 
distribution in the limit of a large number of channels in the lead f42i. which is also 
derived from a cubic equation. Our result in (j70p should then be the Hilbert transform 
of the result in [42] . 

6. Conclusions 

Adding tunnel barriers to the leads of a chaotic cavity, to better approximate realistic 
experimental settings, leads to a wealth of additional correlated trajectory structures 
that contribute in the semiclassical approach. These centre around encounters which 
partially enter the lead and for transmission through the cavity can be incorporated 
into the previous treatment with perfect leads [iTl [21] relatively easily. However the 
increase in complexity meant we could only obtain closed form generating functions 
for the moments of the transmission eigenvalues at leading order if we restricted our 
attention to the case where the tunneling probabilities are equal for each channel in 
each lead as in dTT]) . Neatly, this is a complementary regime to the leading order 
RMT results which allow different tunneling probabilities in each channel as long as 
the leads are identical ^5] . 

At subleading order, the previous treatment 21 could likewise be easily extended 
as in [Appendix A but apart from a computational expansion for the low moments, a 
generating function for all moments could only be obtained in the even more restricted 
case of equal leads and equal tunneling probability for all channels as in (jA.7p . This 
suggests the same case would need to be considered at higher orders in inverse channel 
number, limiting the general semiclassical approach to low orders or low moments. 
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For reflection quantities, even more diagrammatic possibilities occur but of greater 
importance are the new quantum mechanical phases that occur at the tunnel barrier 
[10] . These phases are essential for the unitarity of the semiclassical approach and we 
introduced new combinatorial structures to be able to account for them. With these 
structures, the leading order moments of the reflection eigenvalues in (|43t could be 
obtained. A connection to RMT arises when making one lead perfect while leaving a 
tunnel barrier in the other lead as in section [4.41 since the probability distribution of 
the reflection eigenvalues for such a cavity has recently been derived [41] . Integrating 
this distribution and analysing the asymptotics of a large number of channels should 
lead back to our moment generating function (|47|) . However in the case of one perfect 
and one tunneling lead, from the unitarity of the scattering matrix, we can actually 
obtain this moment generating function directly from the treatment with two perfect 
leads. This may suggest that the RMT probability distribution for two tunneling leads 
could be much more complex. 

Including energy dependence in the new auxiliary structures allowed us to obtain 
the density of states of Andreev billiards in [Appendix B| as well as the moment 
generating function for the Wigner delay times. The full probability density of 
the delay times in known from RMT [42] and it would be interesting to see if 
asymptotic analyses of its integral, like those performed without tunnel barriers 
[551 15U] , would recreate the generating function in ([70]) . Also treated in the asymptotic 
expansions of [33l [39] were the transport moments of quantum dots connected to 
superconductors and their corresponding symmetry classes [48j . Semiclassically the 
conductance through such systems [49l [50] involves different tree species and the 
auxiliary structures we introduced for the moments of the reflection amplitudes could 
be useful for investigating the higher moments of these superconducting symmetry 
classes. 

The approximations used for the semiclassical diagrammatic rules are valid in the 
regime where the time trajectories spend in the cavity is shorter than the Heisenberg 
time. Without tunnel barriers, this translates to iV > 1 which, seeing as there must 
be at least one channel in each lead to allow transport, is automatically the case (the 
approximations actually seem to hold all the way to = 1 for quantities with a 
single lead). Effects due to correlations above the Heisenberg time can then safely be 
neglected. However, with tunnel barriers the condition becomes pN > 1 so that even 
with a large number of channels, a weak enough tunneling probability p can allow 
Heisenberg time effects to appear. How to treat such long time trajectory correlations 
semiclassically is currently unknown, although they have been treated indirectly in the 
related problem of the spectral statistics of closed chaotic systems using pseudo-orbit 
correlations and resummation imposing unitarity [5H 152] . 

For the diagrammatic rules, the dwell time that trajectories typically spend in 
the cavity should also be much larger than the Ehrenfest time. As the Ehrenfest time 
becomes comparable, additional diagrammatic structures arise and the semiclassical 
evaluations become more complicated. This has limited the treatment without tunnel 
barriers to low moments and low orders [53] [54] [55] [56] [57] although the leading 
order contribution for all moments has also been obtained [58]. With tunnel barriers 
this regime has only been investigated for fewer quantities at lower moments and 
order [401 159] , with particularly interesting behaviour arising from the interplay of the 
tunnel barriers and Ehrenfest time effects for diagrams contributing to the universal 
conductance fluctuations [45] . 
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Figure Al. At subleading order, the trajectory quadruplet contributing to 
tiic second moment in (a) involves a loop traversed in opposite directions by a 
trajectory and its partner and so requires time reversal symmetry. The trajectories 
can be redravi?n around the graph in (b) which is a Mobius strip with an encounter 
node where trees connecting to the channels meet. With tunnel barriers we can 
allow this encounter node to touch the incoming lead and have one (c) or both 
(d) of the trajectory stretches going to the incoming channels tunnel directly into 
the lead. The rest of the encounter must be reflected at the tunnel barrier so that 
the remaining trajectory stretches can travel to the other lead. The encounter 
may also move into the outgoing lead as in (e)-(g) where the reflected trajectory 
stretches may now be along the Mobius loop itself and, as in (g), can become 
separated in the graphical representation. 
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Appendix A. Subleading order for transmission 

Diagrams at subleading order, like the example in figure lAlT a). all involve a loop 
traversed in opposite directions and therefore require time reversal symmetry. This 
loop can be represented graphically as a Mobius strip, while the rest of the diagram 
remains as trees which meet the Mobius strip at encounter nodes [21]. The trajectory 
quadruplet in figure lAlT a) then transform to the graph in figure lAlT b). Each l- 
encounter leads to ^ — 1 of each of the leading order trees cj) and (/). These are arranged 
alternately around the encounter with an arbitrary number on each side of the Mobius 
strip. If the numbers on each side are odd or even, we correspondingly describe the 
node as odd or even. One complications is that to close the Mobius loop we require 
an odd number of odd nodes around it. 

With tunnel barriers, we can also move the encounters into the incoming or 
outgoing leads and allow some links from the encounter nodes to tunnel straight into 
the leads while the remaining parts of the diagram are reflected back into the cavity. 
From the diagram in figure lAlT a) or (b), we obtain the additional possibilities in 
figures EHc)-(g). For the semiclassical contribution B of each odd node, for which 
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there are I — 1 ways of having an odd number of trees on each side, we obtain 

oc N 

f = - EE(i-(i-P')Oa"i)rv-^ 

^ 1=2 1=1 

AT l-l 



1=2 i=Ni + l k=l 

^ / / 1 ^ ^ 



, ^^a-i)'^'-'!^ / l0'-^-nl-P^r'=pfr^ (a.i) 

(=2 1=1 /c=l ^ ^ 

where we include a factor of q to later count the number of odd nodes and do not 
include the contribution from the links which form the Mobius strip itself. Again the 
(1 — pY term corresponds to the k = terms of both sums, leading to 

B _ N4>i ^ y. 0(l-K)[rK + 0(l-p,)] 

^ "P»)[?'P» + 0(l-p»)] ,^2) 

jr[ [l-rp,0-00(l-p,)]' ' 
There are / ways of arranging an even number of trees on each side, giving a further 
contribution of 

N 



j^^o ^lyy ■y 
" 9 1-00 ._tr",, 1 



(A.3) 



^ 0(1 "P»)[rp^ +(j){l-p,)] 



t=i l-''p^0-00(l-p,) 

Around the Mobius strip we arrange an arbitrary number of nodes each separated 
by a link, account for the rotational symmetry by dividing by the number of nodes, 
and then ensure we have an odd number of odd nodes [21] to obtain the generating 
function 

ki = -llnil-yiA + B)], ^ K.i^ - D ^ -I) ^ ^^^^^ 

which generates all diagrams without fixing any of the channels to be the first. This 
freedom allows an additional factor of 2n which can be obtained by differentiating to 
give the desired subleading order moment generating function 

(A.5) 

dr 

However, even when all of the tunneling probabilities are the same, where we can 
use ([13]) and (IT4)) to simplify ()A.2p and (jA.Sp , since and are determined by quartic 
equations we were unable to obtain the corresponding quartic for Ti . Of course we 
can substitute the expansions of and from ([19]) into ()A.2)) and ([A.SP and obtain 
arbitrarily many moments in the expansion of Ti : 

Ti = - ^ps + 2^p (1 - 2p - 2^(3 - 4p)) s2 

+ ^p (3p(2 - 3p) + C(18 - 93p + 87p2) _ ^{90 - 279p + 205^^)) + . . . (A.6) 

The first two terms have been obtained previously [301 1131 Si] , for the more general 
case of arbitrary tunneling probabilities which could likewise be incorporated in 
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the expansions here. The first term which is the weak localisation correction to 
the conductance agrees with the known semiclassical [IHl SI] and RMT [35] result 
with equal probabilities, while the correction to the second moment matches the 
corresponding semiclassical result from considering the diagrams explicitly [45] which 
in turn is equal to the RMT diagrammatic result [43) . 

In the further simplified case of equal leads = (^2 where 4> = (j) are given by the 
quadratic in ([22p . we can directly obtain the corresponding generating function 

sp{l~-p){2 - p)(2 - s)VT~s-sp\l -p){2- 2s + s^)~ s^p^A 
' [4(l-p) + sp2][4(l-p)(s-l) + sp2](l-s) ■ 

Appendix B. Andreev billiards 

Andreev billiards are systems where the (single) normal conducting lead has 
been replaced by a superconductor to make a closed ballistic system trapping 
electrons which are converted to holes (and vice versa) when they interact with the 
superconductor. Interference effects, as exemplified by encounters between trajectory 
sets in the semiclassical tree recursions, lead to an energy region where the system 
supports no quantum states and a hard gap in the density of states [THl 110] • The 
semiclassical treatment of Andreev billiards was detailed in [50] so we just highlight 
the incorporation of tunnel barriers here. 

In the scattering approach [61_ , the density of states [S5] 

n oe 

n— 1 

involves a slowly converging series of the correlation functions C(e, n) defined in (IMl) 
and generated in section [SJ For the density of states however we would rather directly 
obtain (again with a — ie) 

-'^ n oa 

n— 1 

where to divide by n we can generate the trees without rooting them in a specific 
channel, as was automatically the case when looking beyond leading order in 
[Appendix A[ At leading order [^ , we can join 21 subtrees at a single point to make 
an encounter and divide by the rotational symmetry factor of 21, 

r> Y^(l-^a)/^' , . If;.! + 2w;.o + i&;.o 

-2. — 21 — + ^2.- ^ 

1=2 1=2 



((1 - p)^wi,i - 2p{l - p)wifl + p^wifi/{l - p)) 

, (B.3) 



to overcount all unrooted trees by the number of their encounter nodes. The last two 
terms correspond to allowing the encounter to touch the lead with the top link either 
being reflected or tunneling straight through, while the (1— p)' term from the encounter 
is still included in w;,o- To perform the sums divided by I we can include powers of 
q in the recursion relations for w and w in psp and find and solve the corresponding 
coupled equations for '}2d^2 (z'^^wj.iand X](^2 'z'^^^i.o- Integrating with respect to q 
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Figure Bl. The leading order contribution to the density of states of Andreev 
billiards with a tunnel barrier between the superconductor and the cavity for 
different tunneling probabilities p. 



and setting it equal to 1 then provides the required sums and we obtain 

2ko = ln(l - f) + - In [1 + (p - l)(r2 + f) + rf{rf - 2p)] 

+ (p-l)r2 + /p(/-2r). (B.4) 
We can also fuse two subtrees (with at least one encounter node each) together, 

= ip-a)(f- ' - - (1 - p)r^ (B.5) 

to overcount all unrooted trees by the number of their internal edges, which is one less 
than the number of encounter nodes. The last two terms are corrections for the first 
moment. The difference Kq — Kq then counts the unrooted trees exactly once and we 
obtain the desired generating function 

which provides the leading order terms of C{e,n)/n. Differentiating implicitly with 
respect to the energy difference, or a, gives the quartic satisfied by Hq 
ap^{s - if {aH^ - 2Hl) + a^pil - s)(2 + p + 2s - isp)Hl 

+ [(1 - s){p + a)+pasf {Hi + iJo) + 2pa{s - l){2 + p - ps - a - as + pas)Hl 

+ 4pa{s - l)Ho - sp^ = 0. (B.7) 

We can now make the substitution 

iW^o(e) ~ 1 

Ho{s = -1) = , (B.8) 

and replace a = ie to obtain the leading order contribution to the density of states 
do{e) = — ImWo(e) through the quartic 

p'^e^W;^ + Ap^eWl^ + {Ap'^+p'^e^ + Ape^ - Ae^) 

+ Ap^eWo + Ape^ - Ae^ = 0. (B.9) 
This is identical to the RMT result |60j and some example density of states for various 
transparencies p of the tunnel barrier are depicted in figure IBll 
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